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ABSTRACT 

Assuming that damped Lyman-a (DLA) systems are galactic discs, we calculate 
the corresponding evolution of metal abundances. We use detailed multi-zone models 
of galactic chemical evolution (reproducing successfully the observed properties of 
disc galaxies) and appropriate statistics (including geometrical propability factors) to 
calculate the average metallicity as a function of rcdshift. The results arc compatible 
with available observations, provided that observational biases are taken into account, 
as suggested by Boisse et al. (1998). In particular, high column density and high 
metallicity systems are not detected because the light of backround quasars is severely 
extinguished, while low column density and low metallicity systems are not detectable 
through their absorption lines by current surveys. We show that these observational 
constraints lead to a "no-evolution" picture for the DLA metallicity, which does not 
allow to draw strong conclusions about the nature of those systems or about their role 
in "cosmic chemical evolution" . 
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1 INTRODUCTION 

Damped Lyman-a systems (DLAs) are high column density 
(N(HI)>10'^'' cm~'^) absorbers detected in the optical spec- 
tra of quasars up to relatively high redshifts (up to z ~5). 
Their study constitutes a powerful means to investigate the 
properties of distant galaxies (or of their building blocks). 
In particular, DLA metal abundances have been widely used 
in the past few years (e.g. Prochaska and Wolfe 1999, Pet- 
tini et al. 1999, Edmunds and Philips 1997, and references 
therein) in order to probe the nature of DLAs (i.e. galac- 
tic discs vs. galaxies of other morphological types) or even 
to probe the so-called "cosmic chemical evolution" (i.e. the 
global evolution of gas, metallicity and star formation rate 
in the Universe). It is not clear, however, whether the ob- 
served abundances alone allow to probe the nature or the 
history of those systems. One reason is possible depletion of 
metals into dust (e.g. Pei and Fall 1995). In this work we in- 
vestigate another factor, namely observational biases, along 
the lines suggested by Boisse et al. (1998). Assuming that 
(proto)galactic discs constitute a quite plausible model for 
such systems we show that observational biases may com- 
pletely alter our interpretation of DLA metal abundances. 



2 OBSERVATIONAL BIASES OF METAL 
ABUNDANCE ESTIMATES IN DLAS 

In Fig. 1 we present the empirical evidence on which this 
work is based. Observations of Zn abundances in DLAs ob- 



tained by several groups are plotted as a function of HI 
column density N(HI). This plot is an up-dated version of 
Fig. 19 of Boisse et al. (1998) and confirms the suggestion 
of those authors, namely that there seems to be an anti- 
correlation between the observed Zn abundance and N(HI), 
independently of redshift z. As Boisse at al. (1998) notice, 
this is not to be interpreted as a physical correlation (i.e 
that high metallicities are characteristic of low column den- 
sity systems); indeed, the observed abundance gradients in 
spiral galaxies offer clear evidence that higher abundances 
are found in the inner disc, where gas column densities are 
also higher (e.g. Garnett et al. 1997). 

The correlation of Fig. 1 is rather to be interpreted in 
terms of observational biases: no systems with a combina- 
tion of metallicity and column density outside the shaded 
region of Fig. 1 are presently detected, even if such sys- 
tems do exist. The lack of high metallicity and high column 
density systems should be attributed to extinction effects, 
since extinction depends on both those factors. The lack of 
low metallicity and low column density systems is attributed 
by Boisse et al. (1998) again to observational selection ef- 
fects, i.e. to the fact that below some level the amount of 
Zn atoms along the line of sight is insufficient to allow for a 
proper spectroscopic detection. 

How "realistic" are these empirically determined de- 
tectability constraints and how can they be quantitatively 
understood? We postpone a discussion of those matters to 
Sec. 4. We simply wish to illustrate here the effect of these 
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constraints (assuming they are real) on the detectability of 
a Milky Way type disc through its metal absorption lines. 

In a recent work (Boissier and Prantzos 1999, herefter 
BP99) we presented a detailed model for the chemical and 
spectrophotometric evolution of the disc of our Galaxy. The 
galactic disc is simulated as an ensemble of concentric, inde- 
pendently evolving rings, gradually built up by infall of pri- 
mordial composition. The chemical evolution of each zone is 
followed by solving the appropriate set of integro-differential 
equations. The spectrophotometric evolution is followed in 
a self-consistent way, i.e. with the star formation rate '^{t) 
and metallicity Z{t) of each zone determined by the chemical 
evolution, and the same stellar Initial Mass Function (from 
Kroupa et al. 1993). The adopted stellar yields, lifetimes, 
evolutionary tracks and spectra are metallicity dependent. 
Dust absorption is included according to the prescriptions 
of Guiderdoni et al. (1998) and assuming a "sandwich" con- 
figuration for the stars and dust layers (Calzetti et al. 1994). 
The star formation rate (SFR) is locally given by a Schmidt- 
type law, i.e proportional to some power of the gas surface 
density Eg and varies with galactocentric radius R as: 

vE'(t,i?) = a.T.g{t,Rf-'' V{R) R-^ (1) 

where V{R) is the circular velocity at radius R. This radial 
dependence of the SFR is suggested by the theory of star 
formation induced by density waves in spiral galaxies (e.g. 
Wyse and Silk 1989). 

It turns out that the number of observables reproduced 
by the model is much larger than the number of free pa- 
rameters. In particular the model reproduces present day 
"global" properties (amount of gas, SFR, and supernova 
rates), as well as the current Milky Way disc luminosities 
in various wavelength bands and the corresponding radial 
profiles of gas, stars, SFR and metal abundances; moreover, 
the adopted inside-out star forming scheme leads to a scale- 
length of ~4 kpc in the B-band and ~2.8 kpc in the K-band, 
in agreement with observations (see BP99). 

At this point it should be noticed that, among all metals 
observed in DLAs through absorption line measurements, 
Zn is usually considered to be the most reliable tracer of 
metallicity, because its abundance is expected to suffer little 
from depletion into dust (e.g. Pettini et al. 1997) However, 
from the theoretical point of view, the nucleosynthesis of 
Zn is not well understood. The only known production site 
is massive stars, and the most detailed models of this site 
are those of Woosley and Weaver (1995, WW95), who give 
the only available models with yields of the various isotopes 
as a function of initial stellar metallicity. The Zn yields of 
these models show an unexplained behaviour: the dominant 
Zn isotope, ^"^Zn, is underproduced, while the yields of the 
next two more important isotopes, ®^Zn and ®*Zn, increase 
strongly as metallicity increases from 0.1 Z© to Z©. The 
reason for this behaviour in the models is not yet well un- 
derstood. As a result, the [Zn/Fe] ratio calculated with the 
metallicity dependent yields of WW95 and a chemical evolu- 
tion model shows a rather abrupt and pronounced increase 
for metallicities [Fe/H]> —1, while it stays nearly constant 
at lower metallicities (Timmes et al. 1995a, Goswami and 
Prantzos 2000). This behaviour is not found in the observed 
pattern of the Zn/Fe ratio: both halo and disk stars of all 
metallicities show solar Zn/Fe (see Goswami and Prantzos 
2000 and references therein) , despite the fact that SNIa, the 
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Figure 1. Zn abundances vs gas column density N(HI). Ob- 
servations in DLAs are from Pettini et al. (1994, 1997, 1999), 
Prochaska and Wolfe (1999), Boisse et al. (1998), Lu et al. (1996). 
Observations are apparently contained within the shaded area, 
which is limited by [Zn/H]+log(N(HI))<21 {upper diagonal) and 
[Zn/H]-|-log(N(HI))>18.8 (lower diagonal). For illustration pur- 
poses, we show the evolution of Zn vs gas column density in three 
different zones of a Milky Way type disc (solid curves, at 0.5, 2.5 
and 5.5 scalelengths from the center, respectively) according to 
our models. The column densities correspond to our model disc 
seen "face-on". The disc evolves from the lower left to the up- 
per right in this diagram, as can be seen from the 3 isochrones 
(dashed curves) at times T=l Gyr, 5.5 Gyr and 13.5 Gyr, respec- 
tively; the corresponding redshifts are for a cosmological model 
with Ho=65 km/s/Mpc, f2o=0.3 and a galaxy formation redshift 
z=6 (i.e. the time of formation of the first stars). 

main Fe producers in the disk, do not produce significant 
Zn amounts. 

This situation does not allow us to trust theoretical nu- 
cleosynthesis prescriptions for the evolution of Zn, since the- 
ory seems unable to reproduce observations in the Galaxy. 
Therefore we adopt an empirical approach for the nucleosyn- 
thesis of Zn, already suggested in Timmes et al. (1995b): 
based on the observed Zn/Fe pattern in the Galaxy, we as- 
sume that Zn traces Fe at all metallicities. Since the history 
of Fe in the Galaxy is observationally well constrained (at 
least in the solar neighborhood) and reasonably well repro- 
duced by our models (see BP99), we assume that it can be 
safely used to trace the history of Zn in other places as well. 

In Fig. 1 we present the evolution of the Zn abundance 
in three representative zones of the Milky Way model (cor- 
responding to the inner disc, the solar neighborhood and the 
outer disc, respectively) as a function of the corresponding 
gas column densities of those zones. It should be noticed 
that in our models we do not distinguish between atomic 
and molecular gas and the curves in Fig. 1 are drawn by 
assuming that all the gas is in atomic form (the contribu- 
tion of He, i.e. 10% by number, is properly removed). This 
approximation certainly affects the results (i.e. the curve 
corresponding to the inner disc, where most of the gas is 
known to be in molecular form, should be shifted to the left 
by 0.5 on that same scale), but for our illustration purposes 
the figure is quite appropriate. 
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Indeed, it is clearly seen that the actual chemical evolu- 
tion of the disc cannot be revealed by observations, because 
of the empirically determined constraints discussed previ- 
ously. At any given epoch (i.e. along the "isochrones" in 
Fig. 1) only a sub-sample of the disc can be probed. At 
early times this sub-sample is representative of the inner re- 
gions which reach rapidely relatively high column densities 
and moderate metallicities; at late times these regions de- 
velop such high densities and metallicities that they move 
into the "forbidden" part of the diagram (upper right). At 
late times then, it is the outer regions that can be probed 
through the absorption lines. 

Obviously, in such conditions it is difficult to derive the 
real chemical history of the system, since only snapshots of 
different regions at different times can be available, not of 
the same regions at different times. 



3 METALLICITY EVOLUTION IN DLAS 

The discussion in the previous section suggests that our pic- 
ture of DLA chemical evolution, obtained through quasar 
absorption line measurements, may be seriously biased. In- 
deed, if taken at face value. Fig. 1 suggests that at early 
times the lowest metallicities and at late times the highest 
metallicities are unobservable. Could this bias modify the 
picture to such an extent as to give the impression of no- 
evolution at all? 

In order to answer quantitatively this question, a work- 
ing model for the evolution of DLAs is needed. Some au- 
thors suggested that DLAs are proto-galactic discs (e.g. 
Prochaska and Wolfe 1996, Ferrini et al. 1997, Prantzos 
and Silk 1998) while others interpreted them as low sur- 
face brightness galaxies (Jimenez et al., 1998), dwarf irreg- 
ulars (Matteucci et al. 1997), galactic halos (Valageas et al. 
1999) or proto-galactic "building blocks" (Haehnelt et al. 
1998, Ledoux et al. 1998). It has also been suggested that 
DLAs differ substantially from the galaxies that contribute 
mostly to the observed SFR in the Universe (e.g. Pettini et 
al. 1999). 

We shall adopt in this work the hypothesis that DLAs 
are galactic discs. We shall use a detailed model we devel- 
oped recently for the chemical and spectrophotometric evo- 
lution of galactic discs (Boissier and Prantzos 2000). It is 
essentially the Milky Way model presented in Sec. 2, ex- 
tended to disc galaxies through "scaling properties" derived 
by Mo, Mao and White (1998) in the framework of the Cold 
Dark Matter (CDM) scenario for galaxy formation. In our 
simplified version of this scenario disc profiles can be ex- 
pressed in terms of only two parameters: maximal rotational 
velocity Vc (measuring the mass of the halo and, by assum- 
ing a constant halo/disc mass ratio, also the mass of the 
disc) and spin parameter A (measuring the specific angular 
momentum of the halo). If all discs are assumed to start 
forming their stars at the same time (but not at the same 
rate!), the profile of a given disc can be expressed in terms 
of the one of our Galaxy. We constructed a grid of 25 mod- 
els caracterised by Vc = 80, 150, 220, 290, 360 km/s and 
X/Xmw = 1/3, 1, 5/3, 7/3, 3, respectively, where Xmw is 
the spin parametre of the Milky Way. Increasing values of 
Vc correspond to more massive discs and increasing values 
of A to more extended ones. The SFR is calculated from Eq. 



1, with the appropriate velocity profile V{R). Notice that 
the efficiency a is not a free parameter, since it is the same 
as in the Milky Way model. 

It turns out that this simple model reproduces fairly 
well most of the main properties of present days discs 
(Boissier and Prantzos 2000): disc sizes and central sur- 
face brigthness, TuUy-Fisher relations in various wavelength 
bands, colour-colour and colour-magnitude relations, gas 
fractions vs. magnitudes and colours, abundances vs. local 
and integrated properties, as well as integrated spectra for 
different galactic rotational velocities. Moreover, as shown 
in Prantzos and Boissier (2000), it also reproduces the ob- 
served abundance gradients in disc galaxies. 

In Fig. 2 we present the results of our models for the Zn 
evolution in DLAs (assumed to be galactic discs) as a func- 
tion of redshift z. For clarity, only the evolution of the inner 
and outer disc is presented in each panel {thin curves, cor- 
responding to zones located at 0.5 and 5.5 scalelengths from 
the centre, respectively). Star formation is assumed to start 
at redshift z=6 for all discs, but any value of z >4-5 would 
produce results similar to those displayed here. The result- 
ing evolution is not very different from that calculated in e.g. 
Prantzos and Silk (1998) for the Milky Way, or in Ferrini et 
al. (1997) with multi-zone disc models. Notice that Malaney 
and Chaboyer (1996), Timmes et al. (1995b), Matteucci et 
al. (1997), Edmunds and Philipps (1997) and Lindner et al. 
(1999) have studied metallicity evolution in DLAs with one- 
zone models, differing by the star formation timescales or 
by the time of the beginning of star formation. 

It can be clearly seen in Fig. 2 that between redshifts 
z ~3 and z ~1 there is substantial metallicity evolution in 
all galactic zones, typically an increase by a factor ~10. Such 
an increase is certainly not observed in the available data, 
which is also displayed on each one of the panels in Fig. 2. 

The region enclosed within thick solid curves in Fig. 2 
is obtained by application of the "empirical constraints" of 
Fig. 1, i.e. by excluding all regions with a combination of 
metallicity and column density F(Zn,N(HI)) = [Zn/H] + 
log(N(HI)) such that F <18.8 or F >21. The resuffing ob- 
servational picture is now completely different from the real 
one: no sizeable evolution in metallicity is observed (except, 
perhaps, in the lowest redshift range, where metallicities are 
somewhat higher than average). 

It should be noticed that the empirical "filter" has been 
applied to our models by assuming that discs are seen "face- 
on" and that all the gas is in the form of HI. These simplify- 
ing assumptions have opposite effects on the derived column 
density along the line of sight: adopting a different inclina- 
tion would increase the column density of our disc models; 
taking into account that part of the gas is in the form of 
H2, would decrease it. The former factor can be treated sta- 
tistically, but not the later. Taking all other uncertainties 
into account (i.e. possible variations in the H2 /HI ratio with 
metallicity, see Combes 1999) we think that Fig. 2 gives a 
rather good first approximation to the real situation. 

The results of Fig. 2 are summarised in the upper panel 
of Fig. 3, wher we plot in the same diagram all the "filtered" 
zones of our models as a function of redshift {shaded area) 
and compare them to observations. The "no-evolution" pic- 
ture is even more clearly seen now, especially in the z ~l-3 
redshift range. A firm prediction of these models is that the 
Zn abundances of DLAs at higher redshifts, in the range 
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Figure 2. Evolution of Zn abundances in our models as a function of redshift z. Results are presented for a grid of 25 disc models, 
caracterised by 5 values of the disc maximal circular velocity Vc (80, 150, 220, 290, 360 km/s, from top to bottom) and 5 values of the 
spin parameter \/\mw (1/3, 1, 5/3, 7/3, 3, from left to right) where Xmw is the corresponding value for the Milky Way disc. For clarity, 
only the evolution of two "extreme" zones is shown for each model, at 0.5 and 5.5 scalelengths from the center {upper and lower curve, 
respectively, in each panel), which span the evolution of the whole disc. The region enclosed within thick curves is obtained by application 
of the empirical selection criterion of Fig. 1 (i.e. 18.8 < [Zn/H]+log(N(HI)) < 21), assuming that all discs are seen "face-on". It can 
be seen that this "filter" selects zones within a restricted range of [Zn/H] value (independently of redshift), and leads to a no-evolution 
picture. Observations (the same in all panels) are from the references listed in the legend of Fig. 1. 



z ~3-5, will be not too different from those already detected 
at z ~l-3. 

Assuming that this first approximation is correct (i.e. 
that DLAS are indeed galactic discs) it is interesting to cal- 
culate the most probable metallicity values expected at a 
given redshift. Even if our basic assumptions are correct, 
this is by no means a trivial task, since it implies the knowl- 
edge of the appropriate statistical factors as a function of 
the redshift. For illustration purposes we adopt here the fol- 
lowing simplified assumptions: 

i) the distribution function of discs in the velocity space 
Vc is time-independent and given by the expression sug- 
gested in Gonzalez et al. 1999 (in the following we simplify, 
for clarity, the notation Vc to V, unless if explicitly stated 
otherwise): 



Fv{V)dV 



v_ 

.V, 



dV 
14 ' 



(2) 



The parameters ^f*, 14, /3 and n are determined in Gon- 
zalez et al. (1999) on the basis of observed TuUy-Fisher re- 
lationships and luminosity (Schechter-type) functions. We 



adopt here the set of parameters of their Table 4 (fifth row, 
LCRS-Courteau data) corresponding to the velocity inter- 
val covered by our models. Our results would not be affected 
much by the choice of another velocity function, since the 
form of Fv always favours discs of low V. 

ii) the distribution function in spin parameter A-space 
is time-independent and given by: 



exp 



ln^{\/X) 



2ol 



dX 

T 



(3) 



with A=0.05 and a\—0.5 (obtained by numerical simula- 
tions, see e.g. Mo, Mao and White, 1998) and Amw=0.06 
for the Milky Way disc (Sommer-Larsen, private communi- 
cation). The A-function favours moderately "compact" discs 
(those with A ~0.04-0.05). 

iii) the probability that a line of sight to a QSO inter- 
cepts a disc in the radius interval [ii, R + dR] is proportional 
to the geometrical cross-section FndR — 2nRdR, favouring 
the detection of the outer regions of the larger discs. 

iv) the distributions Fv, Fx and Fr are independent. 
Applying the joint probability function F{V, A, R) = 
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Figure 3. Evolution of mean metallicities of our disc galaxy 
models as a function of reshift z. Upper panel: The thin curve 
is the mean metallicity of all the zones of our models, i.e. the 
"unfiltered" value < [Zn/H] >u obtained with <I>(R)=1 always 
in Eq. (4); the thick curve is the mean value < [Zn/H] >p of 
the zones "filtered" through the empirical constraint of Fig. 1, 

1. e. the mean value of the shaded area, obtained with ^{R)=l 
for those filtered zones and ^{R)=0 for the others. The shaded 
area includes all the "filtered" zones of our models, i.e. all the 
zones inside thick curves in Fig. 2. Observations are as in Fig. 

2. Lower panel: Column density weighted average metallicity of 
the "filtered zones" < [Zn/H] >FWy obtained by Eq. (5) (thick 
curve); data are from Pettini et al. (1999). 



FvF\Fr to our models, we obtain the mean metallicities 
< [Zn/H] > shown in Fig. 3. The two curves in the upper 
panel are obtained by: 



Zn 



H 



Ix Iv ^(^' V^RMRWiR)] dR dV d\ 



(4) 



where Rl is the radius of the largest disc in our models. 
The mean value over the whole disc (i.e. without apply- 
ing the empirical "filter") corresponds to $(i?)=l in all 
zones and is given by the thin curve in Fig. 3. In that case, 
< [Zn/H] >u {U for "Unfiltered") increases by a factor ~20 
between redshifts 2:=3 and z=l and is clearly below all ob- 
servational data (because the outer, low metallicity, regions 
of the discs are favoured in that case). This shows the im- 
portance of properly taking into account various statistical 
factors, something that has not been done in previous stud- 
ies of DLAs with multi-zone disc models (e.g. in Ferrini et 
al. 1997 and Prantzos and Silk 1998). 

The mean value < [Zn/H] >f over the "filtered' disc 
zones (i.e. those in the shaded aerea of Fig. 3) is obtained 
with "I>(_R)=1 in those zones and <1>(_R)=0 outside them. It is 



shown by the thick curve in Fig. 3. This < [Zn/H] >f value 
is in the lower range of the "filtered" values, again because 
of the geometrical factor Fa. < [Zn/H] >f increases by 
a factor of ~2 between 2:=3 and z=l, an increase which is 
compatible with the observations. 

Pettini et al. (1999) have also estimated the column 
density weighted average of the [Zn/H] values in DLAs, by 
binning their data in 5 redshift bins (the last bin, at z; > 3, 
being in fact an upper limit only). Again, no substantial 
evolution is seen in the data {lower panel in Fig. 3). We 
also calculated the corresponding average metallicity in our 
models, by folding with the gas column densities Nh{R) of 
the "filtered" zones (assuming that all but 10 % by number 
- corresponding to He - is in the form of atomic hydrogen): 



Zn 



H 



Nh dR dV d\ 



/a !v /o ■^(^' ^' '^^ '^'^ 

where $, [Zn/H] and Nh depend on radius R. The resulting 
evolution is shown in the lower panel of Fig. 3. It can be 
seen that the weighted mean metallicity < [Zn/H] >fw oi 
the "filtered" zones of our models evolves in a way which 
is certainly compatible with the data. A clear prediction of 
the model is that at low redshifts {z <1) there should be as 
much evolution as in the 2=1-3 range (i.e. a factor of ~2-3 
increase in the weighted mean metallicity of DLAs in both 
cases) . 

In summary, by using a self-consistent model of galactic 
chemical evolution (i.e. reproducing in detail the properties 
of local galaxies) and incorporating "reasonable" statistics 
and appropriate empirical constraints, we have shown un- 
ambiguously that only a small degree of evolution should be 
expected for the obeserved mean metallicity of DLAs; ac- 
cording to our models, these systems may well be galactic 
discs. 



4 ORIGIN OF BIASES 

It is interesting to see how the empirically determined up- 
per limit in the Zn vs N(HI) plane may be intrepreted in 
terms of extinction. We shall assume here that our model 
discs constitute gaseous screens, seen face-on and reducing 
the intensity of the hght of background quasars. The cor- 
responding extinction is calculated in the rest-frame of the 
absorber by the formula: 

xl.6 



-(5) 



Ax{z) = {Ax/Av)e {Av/Nh)q Nh {Z/ZqY 



(6) 



where {Ax/ Av)q is the local normalised extinction curve 
(Natta and Panagia,1984), {Av/Nh)q is the extinction in 
the V-band in the solar neighborhood (from Bohlin et al., 
1978) and the exponent 1.6 in the metallicity term is intro- 
duced in order to reproduce the observed extinction curves 
in the Small and Large Magellanic Clouds (Guiderdoni and 
Rocca-Volmerange, 1987). The wavelength A depends on the 
redshift z of the absorber. Since the background QSO hght 
is observed in the visible {Xv,obs~0-55 fi), we have for the 
corresponding wavelength Xabs of the absorber 

X^^s{z) = ^ (7) 

For illustration purposes, we calculated the extinction 
at 3 different redshifts 2=1,2,3 (i.e. in the absorber's rest- 
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Figure 4. Upper panel: Extinction of background light 
caused by the various zones of our models, assumed to be gaseous 
screens (seen face-on) at redshift 2=1,2,3 (from left to right, re- 
spectively) of column density N(H) and metallicity [Zn/H] (Eq. 
6) Lower panel: Corresponding reduction factor in the brightness 
of the background light. In both panels the shaded area corre- 
sponds to the one of Fig. 1, i.e. satisfying the "filter" 18.8 < 
[Zn/H]-|-log(N(H)) < 21. It is clearly seen that extinction in- 
creases rapidly to the right of the upper limit, presumably making 
the bakground QSOs unobservable. 



frame Xabs= 0.275, 0.183 and 0.137 fi, respectively) during 
the evolution of all the zones of our models. We plot it in 
Fig. 4 as a function of F=[Zn/H]-f log(N(HI)) (upper panel), 
while in the lower panel we plot the corresponding fraction of 
background light filtered through the screen. It can be seen 
that for values of F >21 (i.e. above the empirically deter- 
mined limit in Fig. 1) extinction increases rapidly, reaching 
1 mag at z=l, 1.5 mag at z=2 and 2 mag at z—3). The 
background intensity drops below 40% of its initial value. 
Our results substantiate the claim of Boisse et al. (1998) 
that extinction is biasing the interpretation of metallicity 
abundance determinations in DLAs. 

As for the lower value of the "empirical" constraint, it 
can be understood as follows: Taking into account that, by 
definition, DLAs correspond to column densities log(N(HI)) 
>20, the lower limit F <18.8 corresponds to [Zn/H]<-1, i.e. 
to less than 5 10^^ atoms of Zn per cm^ along the line of sight 
(adopting a solar ratio (Zn/H)© ~4.5 10~** by number, e.g. 
Anders and Grevesse 1989). Current surveys of DLAs typ- 
ically reach N(Zn)>10^^ cm"^ (Pettini et al. 1997), which 
explains, perhaps, why the lower left part of Fig. 1 is void. 

In any case, it should be interesting to detect the stellar 
population responsible for the chemical enrichment of the 
gas in DLAs. In Fig. 5 we plot our model B-V values of 




3 

Redshift 



Figure 5. B-V vs. redshift for the stellar populations of the 
disc zones of our models that satisfy the condition 18.8 < [Zn/H] 
-f log(N(HI)) < 21 (corresponding to the shaded aerea of Fig. 
3). Discs are assumed to be seen "face-on"; inclined discs should 
provide higher B-V values. 



that population, which span a narrow range 0.2<B-V<0.4 
for redshifts 1< z <3. In fact, these values are lower limits, 
since our model discs are assumed to be seen face-on. Thus, 
one of the firm predictions of our model applied to DLAs is 
that the underlying stellar population should be somewhat 
redder than the curve in Fig. 4. 



5 SUMMARY 

The observationally deter- 

mined constraint 18.8 < [Zn/H]-f log(N(HI)) < 21 (Boisse et 
al. 1998) has profound implications for the intrepretation of 
metal abundances detected in DLAs. Assuming that DLAs 
are galactic discs, and using detailed (and successful) models 
for disc evolution, we show that current observations can- 
not probe the true evolution of those systems: a no-evolution 
picture, compatible with available data, emerges when the 
empirical constraints are taken into account. We calculate 
average metallicities in the galactic zones of our models that 
are "filetered" by the empirical constraints, by taking into 
account appropriate statistical factors. We find that the re- 
sulting column density weighted average metallicity shows a 
small increase at low redshifts and is compatible with cur- 
rently available data. These findings suggest that DLAs may 
well be galactic discs, as argued e.g. in Prochaska and Wolfe, 
1999 (for a different view, see Pettini et al., 1999). 

We also show quantitatively how extinction may be in- 
deed responsible for the non-detection of metal rich DLAs 
of high column density, as suggested by Boisse et al. (1998); 
according to those authors, such DLAs may be one day 
detected in samples drawn from the observation of fainter 
QSOs, independently of the redshift. Moreover, metal poor 
DLAs of low column density should also be detected with 
more sensitive instruments. Finally, if our interpretation of 
currently observed DLAs as galactic discs is correct, we ex- 
pect that the underlying stellar populations should have B- 
V>0.2 (by a small amount) in the redshift range 1< z <3. 
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